Control method of a temperature of a sample

ABSTRACT

A method of stably controlling the temperature of a sample placed on a sample stage to a desired temperature by estimating a sample temperature accurately, the sample stage including a refrigerant flow path to cool the sample stage, a heater to heat the sample stage, and a temperature sensor to measure the temperature of the sample stage. This method comprises the steps of: measuring in advance the variation-with-time of supply electric power to the heater, temperature of the sample, and temperature of the temperature sensor, without plasma processing; approximating the relation among the measured values using a simultaneous linear differential equation; estimating a sample temperature from the variation-with-time of sensor temperature y 1 , heater electric power u 1 , and plasma heat input by means of the Luenberger&#39;s states observer based on the simultaneous linear differential equation used for the approximation; and performing a feedback control of sample temperature using the estimated sample temperature.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to a sample temperature control method, and more specifically relates to a method of temperature control in a plasma processing apparatus.

2. Description of the Related Art

FIG. 18 is a diagram illustrating a conventional method of adjusting sample temperatures. As shown in FIG. 18, a sample is placed on a sample stage 10 cooled by a refrigerant. The sample stage has an embedded heater 13 and is configured to regulate an amount of heat by regulating the electric power to be supplied to the heater 13. Furthermore, the sample stage has an embedded temperature sensor 15 capable of measuring a temperature of the sample stage 10 (refer to Japanese Unexamined Patent Application Publication No. 2005-522051 and its corresponding U.S. Pat. No. 6,921,724).

FIG. 19 is a diagram illustrating a conventional temperature control apparatus. As shown in drawing, the apparatus employs a feedback control method for electric power u₁ to be supplied to the heater so that measurement value y₁ of the temperature sensor becomes a target value r. Also, the PI control is commonly used as a feedback control method. In recent years, a structure wherein the sample stage is divided into three regions: center, edge, and middle, each having a heater and a temperature sensor has been realized.

SUMMARY OF THE INVENTION

In the conventional method described above, it is possible to quickly control the temperature of a sample stage. However, (1) heat transfer between the sample stage and a sample placed thereon is not sufficient. Therefore, it takes a long time for wafer temperature to reach a desired value. Also, (2) as described above, in a case where the sample stage is divided into several regions and wafer temperature is controlled for each region, the temperature gradient of the wafer may decrease due to heat transfer of the wafer itself when sloping wafer temperature while differing set values of sensor temperature. In such a case, wafer temperature distribution will not agree with temperature distribution of the sample stage. In addition, (3) during the actual etching process the wafer is heated by a heat input from a plasma. Due to this, wafer temperature rises above a target temperature.

In the conventional method described above, it is possible to quickly control the temperature of a sample stage. However, (1) heat transfer between the sample stage and a sample placed thereon is not sufficient. Therefore, it takes a long time for wafer temperature to reach a desired value. Also, (2) as described above, in a case where the sample stage is divided into several regions and wafer temperature is controlled for each region, the temperature gradient of the wafer may decrease due to heat transfer of the wafer itself when sloping wafer temperature while differing set values of sensor temperature. In such a case, wafer temperature distribution will not agree with temperature distribution of the sample stage. In addition, (3) during the actual etching process the wafer is heated by a heat input from a plasma. Due to this, wafer temperature rises above a target temperature.

The present invention has been worked out in view of the problems described above, and therefore an object of the invention is to provide a sample temperature control method capable of stably controlling a wafer temperature so as to attain a desired temperature by estimating a sample temperature accurately.

In order to solve the problems described above, the present invention employs the following means.

A sample temperature control method for controlling the temperature of a sample placed on a sample stage that is disposed in a plasma processing chamber for placing the sample, and has a refrigerant flow path to cool the sample stage, a heater to heat the sample stage, and a temperature sensor to measure a temperature of the sample stage, the method comprising the steps of: measuring in advance, with no plasma processing being performed, the variation-with-time of supply electric power to the sample stage, temperature of the sample, and temperature of the temperature sensor; approximating the relation among these measured values using a simultaneous linear differential equation; estimating a sample temperatures from the variation-with-time of sensor temperature, heater electric power, and plasma heat input by means of the Luenberger's states observer based on the linear simultaneous differential equation used for the approximation; and performing a feedback control of sample temperature using the estimated sample temperatures.

Since the present invention is configured as described above, it is possible to accurately estimate a sample temperature and stably control the sample temperature to a desired temperature.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the present invention will be described in detail based on the following figures, wherein:

FIG. 1 is a diagram illustrating a first wafer temperature estimation method;

FIG. 2 is a diagram illustrating a second wafer temperature estimation method;

FIG. 3 is a diagram illustrating a third wafer temperature estimation method;

FIG. 4 is a diagram illustrating a method of finding a plasma heat input under actual process conditions;

FIG. 5 is a diagram illustrating a temperature control apparatus to be used in the present invention;

FIG. 6 is a diagram illustrating a microwave plasma etching apparatus;

FIG. 7 is a diagram illustrating the details of a sample stage;

FIG. 8 is a diagram showing the variation of an estimated value of wafer temperature according to the first estimation method;

FIG. 9 is a diagram showing the variation of a measured value of wafer temperature according to the first estimation method;

FIG. 10 is a diagram showing the variation of a sensor temperature according to the first estimation method.

FIG. 11 is a diagram showing the variation of a measured value of wafer temperature according to the first estimation method;

FIG. 12 is a diagram showing the variation of an estimated value of wafer temperature according to the second estimation method;

FIG. 13 is a diagram showing the variation of a measured value of wafer temperature according to the second estimation method;

FIG. 14 is a diagram showing the variation of an estimated value of wafer temperature according to the third estimation method;

FIG. 15 is a diagram showing the variation of a measured value of wafer temperature according to the third estimation method;

FIG. 16 is a diagram showing the variation of an estimated value with a plasma heat input according to the second estimation method;

FIG. 17 is a diagram showing the variation of a measured value with a plasma heat input according to the second estimation method;

FIG. 18 is a diagram illustrating a conventional method of adjusting sample temperatures; and

FIG. 19 is a diagram illustrating a conventional temperature control apparatus.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

Preferred embodiments of the present invention are described below with reference to the accompanying drawings. FIG. 5 is a diagram illustrating a temperature control apparatus to be used in the present invention. As shown in FIG. 5, a continuously varying wafer temperature x₂ is estimated from the variation-with-time of sensor temperature y₁ and heater electric power u₁, and then the estimated value ({tilde over (x)}₂) (a ˜ put over x is expressed as ({tilde over (x)}) for convenience) are fed back to heater electric power u₁ to control wafer temperature x₂. As a wafer temperature estimation method, the following three types of method are used.

FIG. 1 is a diagram illustrating a first wafer temperature estimation method. In this method, at step S1-1, the variation-with-time of wafer temperature, sensor temperature, and heater electric power is measured by varying the electric power of each heater described below, with a dummy wafer having a temperature measurement mechanism being placed on a sample stage in advance and with plasma being off. At step S1-2, the co-relation among wafer temperature, sensor temperature, and heater electric power is approximated using a linear differential equation. At step S1-3, an estimated value of wafer temperature is calculated from sensor temperature and heater electric power during the processing of an actual sample by means of the Luenberger's states observer using the linear differential equation at step S1-2.

FIG. 2 is a diagram illustrating a second wafer temperature estimation method. First, at step S2-1, the variation-with-time of wafer temperature, sensor temperature, and heater temperature is measured by varying each heater electric power by the same method as in the first method. At step S2-2, a sample stage surface temperature is measured by varying the heater electric power with no sample being placed on the sample stage, in exactly the same way as at step S2-1. At step S2-3, the co-relation among wafer temperature, sensor temperature, heater electric power, and sample stage surface temperature is approximated using a linear differential equation. At step S2-4, estimated values of sample stage surface temperature and wafer temperature are calculated from sensor temperature and heater electric power during the processing of an actual sample by means of the Luenberger's states observer using the linear differential equation at step S2-3.

FIG. 3 is a diagram illustrating a third wafer temperature measurement method. At step S3-1 to 3-4, estimated values of sample stage surface temperature and wafer temperature are calculated by the same method as in the second method. These estimated values are designated as sample stage surface temperature estimated value 1 and wafer temperature estimated value 1 respectively. At step S3-5, the variation-with-time of wafer temperature, sensor temperature and heater electric power is measured by varying a plasma heat input condition and heater electric power, with a dummy wafer being placed on the sample stage and with plasma discharge being on. At step S3-6, the correlation among wafer temperature, sensor temperature, heater electric power, plasma heat input, sample stage surface temperature estimated value 1, and wafer temperature estimated value 1 is approximated using a linear differential equation. At step 3-7, plasma heat input under actual process conditions is determined by the method described below. At step S3-8, sample stage surface temperature estimated value 1 and wafer temperature estimated value 1 as well as an estimated value of actual wafer temperature are calculated from the variation-with-time of sensor temperature and heater electric power as well as plasma heat input by means of the Luenberger's states observer using the linear differential equation at step S3-6.

Now, a method of determining a plasma heat input under actual process conditions is described. FIG. 4 is a diagram illustrating the outline of this method. At step S4-1, wafer temperature and sensor temperature are measured when a steady-state is reached under several plasma heat input conditions, with a dummy sample being placed on the sample stage and with the heater electric power being off. At step S4-2, assuming that the relation between sensor temperature and wafer temperature is linear, a transformation matrix from sensor temperature to wafer temperature is formulated. At step S4-3, a steady-state value of sensor temperature is measured under actual process conditions. At step S4-4, a plasma heat input is calculated from the steady-state value of sensor temperature and the steady-state value of wafer temperature using the linear differential equation at step S306.

Embodiment 1

FIG. 6 is a diagram illustrating a microwave plasma etching apparatus having an automatic wafer temperature control function. In this apparatus, a microwave generated by a magnetron 5 is introduced into a decompression treatment chamber 1 via a wave guide 6 and a quartz window 7 to produce a plasma 8. At this time, a treatment gas introduced from a gas inlet 9 is dissociated by the plasma 8 and radicals and positive ions are produced. Also, a butterfly throttle valve 3 is provided between the decompression treatment chamber 1 and a pump and therefore it is possible to regulate a pressure in the decompression treatment chamber by regulating the degree of opening of the throttle valve 3.

Furthermore, a wafer 11 to be etched is placed on a sample stage 10 connected to a high frequency power supply. By applying an electric power to the sample stage from the high frequency power supply via a matching box, it is possible to generate a negative bias voltage in the wafer. The sample is etched by irradiating to the wafer, ions in the plasma which have been accelerated by a negative bias voltage (hereinafter, this high frequency electric power is referred to as a bias electric power).

Also, a window 12 made of barium fluoride is provided on a side surface of the decompression treatment chamber, and is configured to be capable of measuring a sample stage surface temperature with a radiation thermometer in the absence of a plasma. Furthermore, this apparatus has a circulating refrigerant cooling device configured to cool the entire sample stage by a refrigerant that is cooled to a certain temperature by the circulating refrigerant cooling device and is circulated between the circulating refrigerant cooling device and the sample stage.

FIG. 7 is a diagram illustrating the detail of the sample stage 10. The surface of sample stage 10 is covered with a dielectric film splayed thereon in which a pair of positive and negative ESC electrodes 21 a and 21 b is embedded so that a wafer 11 is attracted to the sample stage 10 by application of a direct current across the electrodes. Also, the dielectric film has embedded heaters 13 a, 13 b, and 13 c disposed respectively in three regions, center, middle, and edge, for controlling the sample stage temperature distribution, and these regions can be heated independently by respective heater power supplies 14 a, 14 b, and 14 c. Furthermore, temperature sensors 15 a, 15 b, and 15 c are embedded in the respective regions of the dielectric film to measure the temperatures of the center, middle, and edge regions of the sample stage, so that outputs of respective heater power supplies are controlled by a processing device 16 based on the output signals of respective temperature sensors. In addition, He with a certain pressure is filled between the rear surface of a wafer and the sample stage in order to improve the heat transfer between the wafer and the sample stage.

With reference to the etching apparatus shown in FIG. 6, a method of performing a first wafer temperature estimation control is described according to the procedure in FIG. 1.

First, at step S1-1, heater electric power u₁ is incrementally increased and then the variation of wafer temperature x₁ and the variation of sensor temperature y₁ are measured. Specifically, refrigerant temperature is set to 30° C. and He pressure is regulated to 1 kPa, with a dummy wafer having the temperature measuring function being placed on the sample stage and attracted thereto. After a sufficient period of time, when wafer temperature has reached to 30° C., outputs u₁ of heater power supplies 14 a, 14 b, and 14 c are increased from 0 W to 1000 W sequentially and the variation-with-time of wafer temperature x₁ and sensor temperature y₁ for each of the center, middle, and edge regions is measured.

At step S1-2, the relation among wafer temperature x₁, sensor temperature y₁, and heater electric power u₁ is approximated using equation (1). Specifically, the values of respective elements of constant matrices A₁₁ to A₂₂, B₁₁ to B₂₁ are determined using the method of least squares.

$\begin{matrix} {{Equation}\mspace{14mu} 1\text{:}} & \; \\ \begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} y_{1} \\ x_{2} \end{bmatrix}} = {{\begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}\begin{bmatrix} {y_{1} - T_{0}} \\ {x_{2} - T_{0}} \end{bmatrix}} + {\begin{bmatrix} B_{11} \\ B_{21} \end{bmatrix}u_{1}}}} \\ {{y_{1} = {\begin{bmatrix} y_{c\; 1} \\ y_{m\; 1} \\ y_{e\; 1} \end{bmatrix}:{{Pt}\mspace{14mu}{sensor}\mspace{14mu}{temperature}}}},} \\ {{x_{2} = {\begin{bmatrix} x_{c\; 2} \\ x_{m\; 2} \\ x_{e\; 2} \end{bmatrix}:{{Wafer}\mspace{14mu}{temperature}}}},} \\ {{{u_{1} = {\begin{bmatrix} u_{c\; 1} \\ u_{m\; 1} \\ u_{e\; 1} \end{bmatrix}:{{Heater}\mspace{14mu}{electric}\mspace{14mu}{power}}}},}\;} \\ {T_{0} = {\begin{bmatrix} t_{0} \\ t_{0} \\ t_{0} \end{bmatrix}:{{Preset}\mspace{14mu}{temperature}\mspace{14mu}{of}\mspace{14mu}{refrigerant}}}} \end{matrix} & (1) \end{matrix}$

Suffixes c, m, and e represent the center, middle, and edge regions respectively.

At step S1-3, wafer temperature is estimated under a desired process condition using the Luenburger's states observer, based on the constant matrix A₁₁ to A₂₂, B₁₁ to B₂₁. Specifically, a variable Z₂ is defined and then the variation-with-time of z₂ is calculated using equation (2).

$\begin{matrix} {{Equation}\mspace{14mu} 2\text{:}} & \; \\ \begin{matrix} {\frac{\mathbb{d}z_{2}}{\mathbb{d}t} = {{P\left( {z_{2} - T_{0}} \right)} + {Q\left( {y_{1} - T_{0}} \right)} + {{PL}\left( {y_{1} - T_{0}} \right)} + {Ru}_{1}}} \\ {{P = \left( {A_{22} - {LA}_{12}} \right)},{Q = \left( {A_{21} - {LA}_{11}} \right)},{R = \left( {B_{21} - {LB}_{11}} \right)}} \\ {{z_{2} = {\begin{bmatrix} z_{c\; 2} \\ z_{m\; 2} \\ z_{e\; 2} \end{bmatrix}:{{States}\mspace{14mu}{observer}}}},{L:{{any}\mspace{14mu}{matrix}\mspace{14mu}{of}\mspace{14mu} 3\mspace{14mu}{rows} \times 3\mspace{14mu}{columns}}}} \end{matrix} & (2) \end{matrix}$

An estimated value ({tilde over (x)}₂) of wafer temperature is determined from the z₂ described above using equation (3).

$\begin{matrix} {{{Equation}\mspace{14mu} 3}:} & \; \\ \begin{matrix} {{\overset{\sim}{x}}_{2} = {z_{2} + {L\left( {y_{1} - T_{0}} \right)}}} \\ {{\overset{\sim}{x}}_{2} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 2} \\ {\overset{\sim}{x}}_{m\; 2} \\ {\overset{\sim}{x}}_{e\; 2} \end{bmatrix}:{{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}}}} \end{matrix} & (3) \end{matrix}$

If calculated with the L set to an appropriate value, the estimated value ({tilde over (x)}₂) will almost agree with wafer temperature when heat input from a plasma is substantially small. Accordingly, if appropriate PI control is performed with respect to the ({tilde over (x)}₂), quick and accurate wafer temperature control is possible. In this embodiment, unit matrix (4) was used as a value of L.

$\begin{matrix} {{Equation}\mspace{14mu} 4} & \; \\ {L = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}} & (4) \end{matrix}$

Using this method, heater electric power was controlled by varying sequentially target values of wafer temperatures of the center/middle/edge regions from 30/30/30° C. to 70/70/70° C. and further to 70/60/50° C. respectively. The variation of the estimated value ({tilde over (x)}₂) of wafer temperature at this time is shown in FIG. 8 and the variation of the measured value x₂ of wafer temperature is shown in FIG. 9. Since wafer temperature varied with the variation of the estimated value ({tilde over (x)}₂), it is proved that wafer temperature can be controlled to a desired temperature quickly.

Next, a conventional method of PI control of heater electric power with respect to sensor temperature was also reviewed for comparison. First, heater electric power u₁ is regulated using a dummy wafer having a temperature measuring function so that the wafer temperatures x₂ of the center/middle/edge regions become (a) 30/30/30° C., (b) 70/70/70° C., and (c) 70/60/50° C., and the sensor temperatures y₁ of the center/middle/edge regions at the time are measured. The measured values were (a) 30/30/30° C., (b) 59/58/55° C., and (c) 59/52/39° C. respectively. These values were set as target values of the sensor temperature of the center/middle/edge regions to control heater electric power. The variation of sensor temperature y₁ at this time is shown in FIG. 10 and the variation of the measured value x₂ is shown in FIG. 11.

Although sensor temperature y₁ has reached the target value in 40 to 70 seconds after starting the process, wafer temperature x₂ has not yet reached the target temperature. That is, it takes a long time to control. Furthermore, in this method, each time process conditions are changed, it is necessary to adjust the heater electric power using a dummy wafer having the temperature measuring function so as to attain a desired wafer temperature and to check the sensor temperature at the time.

As described above, it is proved that wafer temperature can be controlled to a desired value quickly and accurately without having to measure and/or adjust a wafer temperature each time process conditions are changed. Also, even when there is a difference in heat conductivity between the sensor temperature measurement location/sample stage and the sample, quick temperature control without the difference is possible.

Although an example of fixing the value of constant matrix A₁₁ to A₂₂, B₁₁ to B₂₂ to one value was described in this embodiment, in a case where refrigerant temperature is substantially different or refrigerant flow rate is different, an estimated value may differ from a measured value. In such a case, if the matrix under each refrigerant condition is obtained and the matrix is switched simultaneously with the change of a refrigerant, quick and accurate control can be achieved.

Embodiment 2

With reference to the etching apparatus shown in FIG. 6, the method of performing a second wafer temperature estimation control is described according to the procedure in FIG. 2.

First, at step S2-1, the variation-with-time of wafer temperature x₁ and sensor temperature y₁ of each of the center, middle, and edge regions are measured the same as at step S1-1 in the embodiment 1. At step S2-2, heater electric power u₁ is incrementally increased in exactly the same time sequence as at step S2-1, and the variation of sample stage temperature x₃ is measured with a radiation thermometer shown in FIG. 6. At step S2-3, the relation among heater electric power u₁, wafer temperature x₁, sensor temperature y₁, and sample stage temperature x₃, which were measured at steps S2-1 and S2-2, is approximated using equation (5). Specifically, the value of a constant matrix A₁₁ to A₃₃, B₁₁ to B₃₁ is determined using the method of least squares.

$\begin{matrix} {{Equation}\mspace{14mu} 5} & \; \\ \begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} y_{1} \\ x_{2} \\ x_{3} \end{bmatrix}} = {{\begin{bmatrix} A_{11} & A_{12} & A_{13} \\ A_{21} & A_{22} & A_{23} \\ A_{31} & A_{32} & A_{33} \end{bmatrix}\begin{bmatrix} {y_{1} - T_{0}} \\ {x_{2} - T_{0}} \\ {x_{3} - T_{0}} \end{bmatrix}} + {\begin{bmatrix} B_{11} \\ B_{21} \\ B_{31} \end{bmatrix}u_{1}}}} \\ {{y_{1} = {\begin{bmatrix} y_{c\; 1} \\ y_{m\; 1} \\ y_{e\; 1} \end{bmatrix}:{{Pt}\mspace{14mu}{sensor}\mspace{14mu}{temperature}}}},} \\ {{x_{2} = {\begin{bmatrix} x_{c\; 2} \\ x_{m\; 2} \\ x_{e\; 2} \end{bmatrix}:{{Wafer}\mspace{14mu}{temperature}}}},} \\ {x_{3} = {\begin{bmatrix} x_{c\; 3} \\ x_{m\; 3} \\ x_{e\; 3} \end{bmatrix}:{{Sample}\mspace{14mu}{stage}\mspace{14mu}{temperature}}}} \\ {{u_{1} = {\begin{bmatrix} u_{c\; 1} \\ u_{m\; 1} \\ u_{e\; 1} \end{bmatrix}:{{Heater}\mspace{14mu}{electric}\mspace{14mu}{power}}}},} \\ {T_{0} = {\begin{bmatrix} t_{0} \\ t_{0} \\ t_{0} \end{bmatrix}:{{Preset}\mspace{14mu}{temperature}\mspace{14mu}{of}\mspace{14mu}{refrigerant}}}} \end{matrix} & (5) \end{matrix}$

The suffixes c, m, and e represent the center, middle, and edge regions respectively.

Next at step S2-3, wafer temperature under a desired process condition is estimated using the Luenberger's states observer, based on the constant matrix A₁₁ to A₃₃, B₁₁ to B₃₁. Specifically, variables z₂ and z₃ are defined first and then the variation-with-time of z₂ and z₃ is calculated using equation (6).

$\begin{matrix} {{Equation}\mspace{14mu} 6} & \; \\ \begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} z_{2} \\ z_{3} \end{bmatrix}} = {{P\begin{bmatrix} {z_{2} - T_{0}} \\ {z_{3} - T_{0}} \end{bmatrix}} + {Q\left( {y_{1} - T_{0}} \right)} + {P\begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right.} \end{bmatrix}} + {Ru}_{1}}} \\ {P = \begin{bmatrix} {A_{22} - {L_{1}A_{12}}} & {A_{23} - {L_{1}A_{13}}} \\ {A_{32} - {L_{2}A_{12}}} & {A_{33} - {L_{2}A_{13}}} \end{bmatrix}} \\ {{Q = \begin{bmatrix} {A_{21} - {L_{1}A_{11}}} \\ {A_{31} - {L_{2}A_{11}}} \end{bmatrix}},} \\ {R = \begin{bmatrix} {B_{21} - {L_{1}B_{11}}} \\ {B_{31} - {L_{2}B_{11}}} \end{bmatrix}} \\ {{z_{2} = \begin{bmatrix} z_{c\; 2} \\ z_{m\; 2} \\ z_{e\; 2} \end{bmatrix}},} \\ {{z_{3} = {\begin{bmatrix} z_{c\; 3} \\ z_{m\; 3} \\ z_{e\; 3} \end{bmatrix}:{{States}\mspace{14mu}{observer}}}},L_{1},{L_{2}:{{any}\mspace{14mu} 3\mspace{14mu}{rows} \times 3\mspace{14mu}{columns}\mspace{14mu}{matrix}}}} \end{matrix} & (6) \end{matrix}$

The estimated value ({tilde over (x)}₂) of wafer temperature is determined from these z₂ and z₃ using equation (7).

$\begin{matrix} {{Equation}\mspace{14mu} 7} & \; \\ \begin{matrix} {\begin{bmatrix} {\overset{\sim}{x}}_{2} \\ {\overset{\sim}{x}}_{3} \end{bmatrix} = {\begin{bmatrix} z_{2} \\ z_{3} \end{bmatrix} + \begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right)} \end{bmatrix}}} \\ {{{\overset{\sim}{x}}_{2} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 2} \\ {\overset{\sim}{x}}_{m\; 2} \\ {\overset{\sim}{x}}_{e\; 2} \end{bmatrix}:{{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}}}},} \\ {{\overset{\sim}{x}}_{3} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 3} \\ {\overset{\sim}{x}}_{m\; 3} \\ {\overset{\sim}{x}}_{e\; 3} \end{bmatrix}:{{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{sample}\mspace{14mu}{stage}\mspace{14mu}{surface}\mspace{14mu}{temperature}}}} \end{matrix} & (7) \end{matrix}$

If the estimated value is calculated with L₁ and L₂ set to appropriate values, the estimated value ({tilde over (x)}₂) almost agrees with the wafer temperature when heat input from a plasma is sufficiently small. Accordingly, if an appropriate PI control is performed with respect to the estimated value ({tilde over (x)}₂), quick and accurate wafer temperature control is possible. In this embodiment, the unit matrices of equations (8) and (9) were used as the values of L₁ and L₂.

$\begin{matrix} {{Equation}\mspace{14mu} 8} & \; \\ {L_{1} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}} & (8) \\ {{Equation}\mspace{14mu} 9} & \; \\ {L_{2} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}} & (9) \end{matrix}$

Using this method, heater electric power was controlled by changing sequentially the target values of wafer temperature of the center/middle/edge regions from 30/30/30° C. to 70/70/70° C. and further to 70/60/50° C. respectively. The variation of the estimated value ({tilde over (x)}₂) of wafer temperature at the time is shown in FIG. 12 and that of the measured value x₂ of wafer temperature is shown in FIG. 13.

Wafer temperature varies with the variation of the estimated value ({tilde over (x)}₂), which proves that wafer temperature can be controlled to a desired value quickly. Also, it has been proved that wafer temperature can be controlled more quickly and accurately than in the first method shown in FIG. 9.

As described above, using the second wafer temperature estimation control method will allow quicker and more accurate control of wafer temperature as compared with the first wafer temperature estimation control method or any conventional method. Also, in this method, even when there is a difference in heat conductivity between the sensor temperature measurement location/the sample stage and a sample, quick temperature control without the difference is possible.

Although this embodiment was described using an example of fixing the value of matrix A₁₁ to A₃₃, B₁₁ to B₃₁ to one value, in a case where refrigerant temperature is substantially different or refrigerant flow rate is different, an estimated value may differ from a measured value. In such a case, if the matrix under each refrigerant condition is obtained and then the matrix is switched simultaneously with the change of a refrigerant, quick and accurate control can be achieved. Further, although sample stage surface temperature was measured in a vacuum in this embodiment, sufficient results can be obtained even from a measurement in the atmosphere.

Embodiment 3

With reference to the etching apparatus shown in FIG. 6, a third wafer temperature estimation control method is described according to the procedure in FIG. 3. In FIG. 6, at steps S3-1 to S3-4, the same process as in the embodiment 2 is performed and the estimated value ({tilde over (x)}₂) of wafer temperature and the estimated value ({tilde over (x)}₃) of electrode surface temperature are calculated. At step S3-5, the variation-with-time of wafer temperature x₂ and sensor temperature y₁ measured under arbitrary three kinds of plasma conditions and under different heater heating conditions. Specifically, when wafer temperature reached 30° C. after sufficient time has passed since the He pressure was regulated to 1 kPa with a wafer having the temperature measurement function being placed and attracted to the sample stage, the wafer is processed, while controlling the wafer temperature using the method of embodiment 1, under three kinds of conditions: (1) a condition under which the plasma heat input in the center region is large and a condition under which the preset temperatures of the center/middle/edge regions are controlled to 40/30/30° C., (2) a condition under which the heat input in the middle region is large and a condition under which the preset temperatures of the center/middle/edge regions are controlled to 30/40/30° C., and (3) a condition under which the heat input in the edge region is large and a condition under which the preset temperatures of the center/middle/edge regions are controlled to 30/30/40° C., and the variation of wafer temperature x₂ and sensor temperature y₁ is measured.

Next, at step S-6, the relation among wafer temperature x₂, sensor temperature y₁, the estimated value ({tilde over (x)}₂) of wafer temperature calculated using equation (7), the estimated value ({tilde over (x)}₃) of electrode temperature, and plasma heat input u₂ is approximated using equation (10). Specifically, a constant matrix A₄₁ to A₄₄, B₄₁ is determined using the method of least squares.

$\begin{matrix} {{Equation}\mspace{14mu} 10} & \; \\ \begin{matrix} {\frac{\mathbb{d}x_{2}}{\mathbb{d}t} = {{A_{41}\left( {y_{1} - T_{0}} \right)} + {A_{42}\left( {{\overset{\sim}{x}}_{2} - T_{0}} \right)} +}} \\ {{A_{43}\left( {{\overset{\sim}{x}}_{3} - T_{0}} \right)} + {A_{44}\left( {x_{2} - T_{0}} \right)} + {B_{42}u_{2}}} \\ {u_{2} = {\begin{bmatrix} u_{c\; 2} \\ u_{m\; 2} \\ u_{e\; 2} \end{bmatrix}:{{Amount}\mspace{14mu}{of}\mspace{14mu}{plasma}\mspace{14mu}{heat}\mspace{14mu}{input}}}} \end{matrix} & (10) \end{matrix}$

At this time, it is not necessary to accurately measure the amount of plasma heat input u₂ under each condition, and any three types of independent matrices may be used as shown in equation (1).

$\begin{matrix} {{Equation}\mspace{14mu} 11} & \; \\ \begin{matrix} {{{u_{2}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu} 1} = \begin{bmatrix} 1 \\ 0 \\ 0 \end{bmatrix}},} \\ {{{u_{2}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu} 2} = \begin{bmatrix} 0 \\ 1 \\ 0 \end{bmatrix}},} \\ {{u_{2}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu} 3} = \begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix}} \end{matrix} & (11) \end{matrix}$

Next, after the amount of plasma heat input u₂ under an actual process condition is calculated in the way described below at step S3-7, the estimated value ({circumflex over (x)}₂) of wafer temperature in plasma heat input is measured at step S3-8 using the Luenberger's states observer (a ^ put over x is represented as ({circumflex over (x)}) for convenience). Specifically, the estimated value ({circumflex over (x)}₂) of wafer temperature in plasma heat input is calculated using equations (12) and (13). However, the value of the amount of plasma heat input u₂ is changed from 0 to a value obtained at step S3-7 simultaneously at the start of plasma heat input.

$\begin{matrix} {{Equation}\mspace{14mu} 12} & \; \\ \begin{matrix} {\frac{\mathbb{d}z_{4}}{\mathbb{d}t} = {{P^{\prime}\begin{bmatrix} {z_{2} - T_{0}} \\ {z_{3} - T_{0}} \\ {z_{4} - T_{0}} \end{bmatrix}} + {Q^{\prime}\left( {y_{1} - T_{0}} \right)} + {P^{\prime}\begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right)} \\ {L_{3}\left( {y_{1} - T_{0}} \right)} \end{bmatrix}} + {R^{\prime}\begin{bmatrix} u_{1} \\ u_{2} \end{bmatrix}}}} \\ {{P^{\prime} = \left\lbrack {A_{42} - {L_{3}A_{12}\mspace{14mu} A_{43}} - {L_{3}A_{13}\mspace{14mu} A_{44}}} \right\rbrack},} \\ {{Q^{\prime} = {A_{41} - {L_{3}A_{11}}}},{R^{\prime} = \left\lbrack {{- L_{3}}B_{11}\mspace{11mu} B_{42}} \right\rbrack}} \\ {{z_{4} = {\begin{bmatrix} z_{c\; 4} \\ z_{m\; 4} \\ z_{e\; 4} \end{bmatrix}:{{States}\mspace{14mu}{observer}}}},{L_{3}:{{any}\mspace{14mu}{matrix}\mspace{14mu}{of}\mspace{14mu} 3\mspace{14mu}{rows} \times 3\mspace{14mu}{columns}}}} \end{matrix} & (12) \\ {{Equation}\mspace{14mu} 13} & \; \\ \begin{matrix} {{\hat{x}}_{2} = {z_{4} + {L_{3}y_{1}}}} \\ {{{\hat{x}}_{2} = {\begin{bmatrix} {\hat{x}}_{c\; 2} \\ {\hat{x}}_{m\; 2} \\ {\hat{x}}_{e\; 2} \end{bmatrix}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}{in}}}\mspace{31mu}} \\ {{~~~~~~~~~~~}{{the}\mspace{14mu}{presence}\mspace{14mu}{of}\mspace{14mu}{plasma}\mspace{14mu}{heat}\mspace{14mu}{input}}} \end{matrix} & (13) \end{matrix}$

In this embodiment, the unit matrix of equation (14) was used as a value of L₃.

$\begin{matrix} {{Equation}\mspace{14mu} 14} & \; \\ {L_{3} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}} & (14) \end{matrix}$

Next, the detail of the method of calculating the amount of plasma heat input u₂ under an actual process condition as shown at step S3-7 is described according to the procedure of FIG. 4.

At step S4-1, respective processes under the conditions (1), (2), and (3) described above are performed without heater control, wafer temperature and sensor temperature are measured when steady-state is reached, and the difference (amount of change) from the temperature before plasma heat input is calculated. Then, at step S4-2, a transformation matrix from sensor temperature to wafer temperature is determined by means of equation (15) using the amount of change in wafer temperature and sensor temperature.

$\begin{matrix} {{Equation}\mspace{14mu} 15} & \; \\ {{C = {{{\begin{bmatrix} {\Delta\; x_{c\; 1}} & {\Delta\; x_{c\; 2}} & {\Delta\; x_{c\; 3}} \\ {\Delta\; x_{m\; 1}} & {\Delta\; x_{m\; 2}} & {\Delta\; x_{m\; 3}} \\ {\Delta\; x_{e\; 1}} & {\Delta\; x_{e\; 2}} & {\Delta\; x_{e\; 3}} \end{bmatrix}\begin{bmatrix} {\Delta\; y_{c\; 1}} & {\Delta\; y_{c\; 2}} & {\Delta\; y_{c\; 3}} \\ {\Delta\; y_{m\; 1}} & {\Delta\; y_{m\; 2}} & {\Delta\; y_{m\; 3}} \\ {\Delta\; y_{e\; 1}} & {\Delta\; y_{e\; 2}} & {\Delta\; y_{e\; 3}} \end{bmatrix}}^{- 1}\begin{bmatrix} {\Delta\; x_{c\; 1}} \\ {\Delta\; x_{m\; 1}} \\ {\Delta\; x_{e\; 1}} \end{bmatrix}}\mspace{14mu}{{and}\mspace{14mu}\begin{bmatrix} {\Delta\; y_{c\; 1}} \\ {\Delta\; y_{m\; 1}} \\ {\Delta\; y_{e\; 1}} \end{bmatrix}}\mspace{14mu}{are}\mspace{14mu}{amounts}\mspace{14mu}{of}}}{{change}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu}(1)\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}{and}}\mspace{14mu}{{sensor}\mspace{14mu}{temperature}\mspace{14mu}{{respectively}.\begin{bmatrix} {\Delta\; x_{c\; 2}} \\ {\Delta\; x_{m\; 2}} \\ {\Delta\; x_{e\; 2}} \end{bmatrix}}\mspace{14mu}{{and}\mspace{14mu}\begin{bmatrix} {\Delta\; y_{c\; 2}} \\ {\Delta\; y_{m\; 2}} \\ {\Delta\; y_{e\; 2}} \end{bmatrix}}\mspace{14mu}{are}\mspace{14mu}{amounts}\mspace{14mu}{of}}{{change}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu}(2)\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}{and}}{{sensor}\mspace{14mu}{temperature}\mspace{14mu}{{respectively}.\begin{bmatrix} {\Delta\; x_{c\; 3}} \\ {\Delta\; x_{m\; 3}} \\ {\Delta\; x_{e\; 3}} \end{bmatrix}}\mspace{14mu}{{and}\mspace{14mu}\begin{bmatrix} {\Delta\; y_{c\; 3}} \\ {\Delta\; y_{m\; 3}} \\ {\Delta\; y_{e\; 3}} \end{bmatrix}}\mspace{14mu}{are}\mspace{14mu}{amounts}\mspace{14mu}{of}}{{change}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu}(3)\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}{and}}{{sensor}\mspace{14mu}{temperature}\mspace{14mu}{{respectively}.}}} & (15) \end{matrix}$

Next, at step S4-3, a plasma processing is performed under a desired process condition without heater control, and the difference Δy between sensor temperature and initial temperature is determined at the time when steady-state is reached. Then, at step S4-4, the temperature increase Δx caused by plasma heat input is calculated from the temperature difference Δy using equation (16).

$\begin{matrix} {{Equation}\mspace{14mu} 16} & \; \\ \begin{matrix} \begin{matrix} {{\Delta\; x} = {C\;\Delta\; y}} \\ {{{{Increase}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}\Delta\; x} = \begin{bmatrix} {\Delta\; x_{c}} \\ {\Delta\; x_{m}} \\ {\Delta\; x_{e}} \end{bmatrix}},} \end{matrix} \\ {{{Increase}\mspace{14mu}{of}\mspace{14mu}{sensor}\mspace{14mu}{temperature}\mspace{14mu}\Delta\; y} = \begin{bmatrix} {\Delta\; y_{c}} \\ {\Delta\; y_{m}} \\ {\Delta\; y_{e}} \end{bmatrix}} \end{matrix} & (16) \end{matrix}$

Next, at step S4-5, the increase Δy of sensor temperature and the temperature increase Δx due to plasma heat input are substituted in equation (10), and the amount of plasma heat input u₂ under a desired plasma condition is determined assuming that the left side differential term is 0.

Since wafer temperature can be accurately estimated by this method even in the presence of plasma heat input, it is possible to control wafer temperature quickly and accurately by performing an appropriate PI control on the estimated value ({circumflex over (x)}₂).

Using this method, the target temperatures of the center/middle/edge regions were controlled from 30/30/30° C. to 70/70/70° C. and further to 70/60/50° C. The variation of the estimated value ({circumflex over (x)}₂) of wafer temperature at the time is shown in FIG. 14 and that of the measured value x₂ of wafer temperature is shown in FIG. 15. For comparison, the variation of the estimated value ({tilde over (x)}₂) and that of the measured value x₂ when the control is performed using the second wafer temperature estimation control method of the present invention are shown in FIGS. 16 and 17 respectively.

In the second wafer temperature estimation control method, since the measured value x₂ increases above the estimated value ({tilde over (x)}₂) in plasma heat input, 50 to 70 seconds and 110 to 130 seconds after the start of the process, wafer temperature rises above the target temperature. In the third wafer temperature estimation control method, however, the measured value x₂ agrees with the estimated value ({circumflex over (x)}₂) even in the presence of plasma heat input, and therefore it is proved that wafer temperature x₂ can be controlled quickly and accurately by controlling the estimated value ({circumflex over (x)}₂). Also in this method, even when there is a difference in heat conductivity between the sensor temperature measurement location/the sample stage and the sample, a quick temperature control without the difference is possible.

Although an example of fixing the values of a matrix A₄₁ to A₄₄ and B₄₁ to one value was described in this embodiment, there is a case where the estimated value differs from the measured value when the temperature of a refrigerant is substantially different or the flow rate of a refrigerant is different. In such a case, if the matrix is obtained under each refrigerant condition and the matrix is switched simultaneously with the change of a refrigerant, quick and accurate control can be achieved. Furthermore, although temperature of the sample stage was measured in a vacuum in this embodiment, substantial results can be obtained even if measurement is made in the atmosphere.

As described above, according to the embodiments of the present invention, measurements of heater electric power, wafer temperature, and sensor temperature are performed; the relation among them is approximated using a simultaneous linear differential equation; wafer temperature is estimated by means of the Luenberger's states observer using the simultaneous linear differential equation; and a feedback control is performed using the estimated wafer temperature. This allows quick and stable control of wafer temperature. Also, according to the first wafer temperature estimation method, wafer temperature can be estimated accurately when plasma heat input is substantially small. Further, according to the second wafer temperature estimation method, wafer temperature can be estimated more accurately than in the first method. Further, according to the third wafer temperature estimation method, wafer temperature can be estimated accurately even in the presence of plasma heat input. If the feedback control of heater electric power is performed based on the estimated value of wafer temperature obtained in this way, it is possible to control the wafer temperature quickly and stably as described above. 

1. A control method for controlling the temperature of a sample placed on a sample stage disposed in a plasma processing chamber, the sample stage including a refrigerant flow path to cool the sample stage, a heater to heat the sample stage, and a temperature sensor to measure a temperature of the sample stage, the control method comprising the steps of: measuring, in advance of plasma processing being performed, variations-with-time of: heater electric power, temperature of the sample, and temperature of the temperature sensor; approximating a relationship between the measured values by using a simultaneous linear differential equation; estimating a temperature of the sample from the variation-with-time of: sensor temperature, heater electric power, and plasma heat input, by means of Luenberger states observer equations based on at least one value produced by the simultaneous linear differential equation used for the approximation; and performing a feedback control of temperature of the sample, by using the estimated temperature of the sample to perform a PI control.
 2. The control method according to claim 1, further comprising the steps of: measuring, when the temperature of the refrigerant changes value or when the flow rate of the refrigerant changes value, variation-with-time of: heater electric power, temperature of the sample, and temperature of the temperature sensor; formulating a simultaneous linear differential equation after each change in value; and switching the simultaneous linear differential equation to be used for the Luenberger states observer equations, after each change in value of the refrigerant.
 3. The control method according to claim 1, further comprising the steps of: measuring, in the presence of plasma heat input and in the absence thereof, variation-with-time of: heater electric power, temperature of the sample, and temperature of the temperature sensor; estimating a temperature of the sample by using a relationship between them; and performing a feedback control of temperature of the sample, by using the estimated temperature of the sample to perform a PI control.
 4. The control method of claim 1, wherein the simultaneous linear differential equation includes: ${{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} y_{1} \\ x_{2} \end{bmatrix}} = {{\begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}\begin{bmatrix} {y_{1} - T_{0}} \\ {x_{2} - T_{0}} \end{bmatrix}} + {\begin{bmatrix} B_{11} \\ B_{21} \end{bmatrix}u_{1}}}};$ wherein x₂, y₁, u₁, and T₀ are defined as follows, and wherein suffixes c, m, and e represent the center, middle, and edge regions respectively: ${y_{1} = {\begin{bmatrix} y_{c\; 1} \\ y_{m\; 1} \\ y_{e\; 1} \end{bmatrix}:{{Sensor}\mspace{14mu}{temperature}}}},{x_{2} = {\begin{bmatrix} x_{c\; 2} \\ x_{m\; 2} \\ x_{e\; 2} \end{bmatrix}:{{Temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}}}},{u_{1} = {\begin{bmatrix} u_{c\; 1} \\ u_{m\; 1} \\ u_{e\; 1} \end{bmatrix}:{{Heater}\mspace{14mu}{electric}\mspace{14mu}{power}}}},{{T_{0} = {\begin{bmatrix} t_{0} \\ t_{0} \\ t_{0} \end{bmatrix}:{{Preset}\mspace{14mu}{temperature}\mspace{14mu}{of}\mspace{14mu}{refrigerant}}}};}$ and wherein the values of respective elements of constant matrices A₁₁ to A₂₂ and B₁₁ to B₂₂ are determined using the method of least squares.
 5. The control method of claim 4, wherein the Luenberger states observer equations include: $\frac{\mathbb{d}z_{2}}{\mathbb{d}t} = {{P\left( {z_{2} - T_{0}} \right)} + {Q\left( {y_{1} - T_{0}} \right)} + {{PL}\left( {y_{1} - T_{0}} \right)} + {{Ru}_{1}\mspace{14mu}{and}}}$ P = (A₂₂ − LA₁₂), Q = (A₂₁ − LA₁₁), R = (B₂₁ − LB₁₁);  wherein ${z_{2} = {\begin{bmatrix} z_{c\; 2} \\ z_{m\; 2} \\ z_{e\; 2} \end{bmatrix}:{{States}\mspace{14mu}{observer}}}},{{{and}\mspace{14mu}{L:{{any}\mspace{14mu}{matrix}\mspace{14mu}{of}\mspace{14mu} 3\mspace{14mu}{rows}\mspace{11mu} \times \; 3\mspace{14mu}{columns}}}};}$ wherein from the values of variable z₂ that are calculated in the Luenberger states observer equations, an estimated temperature of the sample ({tilde over (x)}₂) is determined with the following equations: ${{\overset{\sim}{x}}_{2} = {z_{2} + {L\left( {y_{1} - T_{0}} \right)}}},{and}$ ${{\overset{\sim}{x}}_{2} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 2} \\ {\overset{\sim}{x}}_{m\; 2} \\ {\overset{\sim}{x}}_{e\; 2} \end{bmatrix}:{{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}}}};$ ${{wherein}\mspace{14mu}{the}\mspace{14mu}{following}\mspace{14mu}{unit}\mspace{14mu}{matrix}\mspace{14mu}{is}\mspace{14mu}{used}\mspace{14mu}{as}\mspace{14mu}{the}\mspace{14mu}{valve}\mspace{14mu}{of}\mspace{14mu}{L:L}} = {\begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}.}$
 6. A control method for controlling temperature of a sample placed on a sample stage disposed in a plasma processing chamber, the sample stage including a refrigerant flow path to cool the sample stage, a heater to heat the sample stage, and a temperature sensor to measure a temperature of the sample stage, the control method comprising the steps of: measuring, in advance of plasma processing being performed, variation-with-time of: heater electric power, temperature of the sample, temperature of the temperature sensor, and temperature of a surface of the sample stage; approximating a relationship between the measured values by using a simultaneous linear differential equation; estimating a temperature of the sample from the variation-with-time of: sensor temperature, heater electric power, and plasma heat input, by means of Luenberger states observer equations based on at least one value produced by the simultaneous linear differential equation used for the approximation; and performing a feedback control of temperature of the sample, by using the estimated temperature of the sample to perform a PI control.
 7. The control method according to claim 6, further comprising the steps of: measuring, when the temperature of the refrigerant changes value or when the flow rate of the refrigerant changes value, variation-with-time of: heater electric power, temperature of the sample, and temperature of the temperature sensor; formulating a simultaneous linear differential equation after each change in value; and switching the simultaneous linear differential equation to be used for Luenberger states observer equations, after each change in value of temperature or flow rate of the refrigerant.
 8. The control method of claim 6, wherein variation of a sample stage temperature is measured; wherein the simultaneous linear differential equation includes: ${{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} y_{1} \\ x_{2} \\ x_{3} \end{bmatrix}} = {{\begin{bmatrix} A_{11} & A_{12} & A_{13} \\ A_{21} & A_{22} & A_{23} \\ A_{31} & A_{32} & A_{33} \end{bmatrix}\begin{bmatrix} {y_{1} - T_{0}} \\ {x_{2} - T_{0}} \\ {x_{3} - T_{0}} \end{bmatrix}} + {\begin{bmatrix} B_{11} \\ B_{21} \\ B_{31} \end{bmatrix}u_{1}}}},{wherein}$ ${y_{1} = {\begin{bmatrix} y_{c\; 1} \\ y_{m\; 1} \\ y_{e\; 1} \end{bmatrix}:{{{Se}{nsor}}\mspace{14mu}{temperature}}}},{x_{2} = {\begin{bmatrix} x_{c\; 2} \\ x_{m\; 2} \\ x_{e\; 2} \end{bmatrix}:{{Temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}}}},\text{}{x_{3} = {\begin{bmatrix} x_{c\; 3} \\ x_{m\; 3} \\ x_{e\; 3} \end{bmatrix}:{{Sample}\mspace{14mu}{stage}\mspace{14mu}{temperature}}}},{u_{1} = {\begin{bmatrix} u_{c\; 1} \\ u_{m\; 1} \\ u_{e\; 1} \end{bmatrix}:\;{{Heater}\mspace{14mu}{electric}\mspace{14mu}{power}}}},{and}$ $\;{{T_{0} = {\begin{bmatrix} t_{0} \\ t_{0} \\ t_{0} \end{bmatrix}:\;{{Preset}\mspace{14mu}{temperature}\mspace{14mu}{of}\mspace{14mu}{refrigerant}}}};}$ wherein suffixes c, m, and e represent the center, middle, and edge regions respectively; and wherein the value of a constant matrix A₁₁ to A₃₃ and B₁₁ to B₃₁ is determined using the method of least squares.
 9. The control method of claim 8, wherein the Luenberger states observer equations include: ${{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} z_{2} \\ z_{3} \end{bmatrix}} = {{P\begin{bmatrix} {z_{2} - T_{0}} \\ {z_{3} - T_{0}} \end{bmatrix}} + {Q\left( {y_{1} - T_{0}} \right)} + {P\begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right)} \end{bmatrix}} + {Ru}_{1}}},{wherein}$ ${P = \begin{bmatrix} {A_{22} - {L_{1}A_{12}}} & {A_{23} - {L_{1}A_{13}}} \\ {A_{32} - {L_{2}A_{12}}} & {A_{33} - {L_{2}A_{13}}} \end{bmatrix}},{Q = \begin{bmatrix} {A_{21} - {L_{1}A_{11}}} \\ {A_{31} - {L_{2}A_{11}}} \end{bmatrix}},{R = \begin{bmatrix} {B_{21} - {L_{1}B_{11}}} \\ {B_{31} - {L_{2}B_{11}}} \end{bmatrix}},{z_{2} = \begin{bmatrix} z_{c\; 2} \\ z_{m\; 2} \\ z_{e\; 2} \end{bmatrix}},{{{and}\mspace{14mu} z_{3}} = {\begin{bmatrix} z_{c\; 3} \\ z_{m\; 3} \\ z_{e\; 3} \end{bmatrix}\text{:}\mspace{14mu}{States}\mspace{14mu}{observer}}},L_{1},{{L_{2}\text{:}\mspace{14mu}{any}\mspace{14mu} 3\mspace{14mu}{rows}\mspace{11mu} \times \mspace{11mu} 3\mspace{14mu}{columns}\mspace{14mu}{matrix}};}$ wherein from the values of variables z₂ and z₃ that are calculated in the Luenberger states observer equations, an estimated temperature of the sample ({tilde over (x)}₂) is determined using the following equation: ${\begin{bmatrix} {\overset{\sim}{x}}_{2} \\ {\overset{\sim}{x}}_{3} \end{bmatrix} = {\begin{bmatrix} z_{2} \\ z_{3} \end{bmatrix} + \begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right)} \end{bmatrix}}},{wherein}$ ${{{\overset{\sim}{x}}_{2} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 2} \\ {\overset{\sim}{x}}_{m\; 2} \\ {\overset{\sim}{x}}_{e\; 2} \end{bmatrix}\text{:}\mspace{14mu}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}}},{and}}\mspace{11mu}$ $\;{{{\overset{\sim}{x}}_{3} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 3} \\ {\overset{\sim}{x}}_{m\; 3} \\ {\overset{\sim}{x}}_{e\; 3} \end{bmatrix}\text{:}\mspace{14mu}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{sample}\mspace{14mu}{stage}\mspace{14mu}{surface}\mspace{14mu}{temperature}}};{and}}\mspace{14mu}$ wherein  the  following  unit matrices  are  used  as  the  value  of  L₁  and  L₂: ${L_{1} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}},{{{and}\mspace{14mu} L_{2}} = {\begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}.}}$
 10. A control method for controlling temperature of a sample placed on a sample stage disposed in a plasma processing chamber, the sample stage including a refrigerant flow path to cool the sample stage, a heater to heat the sample stage, and a temperature sensor to measure temperature of the sample stage, the control method comprising the steps of: measuring, in the presence of plasma heat input and in the absence thereof, variation-with-time of: heater electric power, temperature of the sample, and temperature of the temperature sensor; approximating, for measurement results in the absence of plasma heat input, a relationship between heater electric power, temperature of the sample, and temperature of the sensor by using a first simultaneous linear differential equation; calculating a first estimated value of the temperature of the sample, by means of Luenberger states observer equations based on at least one value produced by the first simultaneous linear differential equation; approximating, for measurement results both in the presence of plasma heat input and in the absence thereof, a relationship between the calculated first estimated value, sensor temperature, temperature of the sample, heater electric power, and plasma heat input by using a second simultaneous linear differential equation; estimating a second estimated value temperature of the sample, from the variation-with-time of the first estimated value, sensor temperature, heater electric power, and plasma heat input, by means of Luenberger states observer equations based on at least one value produced by the second simultaneous linear differential equation; and performing a feedback control by using the estimated temperature of the sample.
 11. The control method according to claim 10, further comprising the steps of: performing a plasma processing; measuring a steady-state value of temperature of the sample; and substituting the steady-state value in the second differential equation, while assuming that the differential term is zero, to calculate a plasma heat input.
 12. The control method according to claim 10, further comprising the steps of: performing a plasma processing; measuring a relationship between steady-state values of sensor temperature and temperature of the sample, assuming that a heater electric power is zero; approximating a first relationship between steady-state values of sensor temperature and temperature of the sample by using the first simultaneous linear differential equation; performing a desired plasma processing; measuring a second steady-state value of sensor temperature, assuming a heater electric power is zero; estimating a second steady-state value of temperature of the sample, from the first simultaneous linear differential equation and the second steady-state value of sensor temperature, for the desired plasma processing; and substituting the second steady-state value in the second differential equation assuming that the differential term thereof is zero, to calculate a plasma heat input.
 13. The control method of claim 6, wherein the first simultaneous linear differential equation includes: ${{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} y_{1} \\ x_{2} \end{bmatrix}} = {{\begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}\begin{bmatrix} {y_{1} - T_{0}} \\ {x_{2} - T_{0}} \end{bmatrix}} + {\begin{bmatrix} B_{11} \\ B_{21} \end{bmatrix}u_{1}}}};$ wherein x₂, y₁, u₁, and T₀ are defined as follows, and wherein suffixes c, m, and e represent the center, middle, and edge regions respectively: ${y_{1} = {\begin{bmatrix} y_{c\; 1} \\ y_{m\; 1} \\ y_{e\; 1} \end{bmatrix}\text{:}\mspace{14mu}{Sensor}\mspace{14mu}{temperature}}},{x_{2} = {\begin{bmatrix} x_{c\; 2} \\ x_{m\; 2} \\ x_{e\; 2} \end{bmatrix}\text{:}\mspace{14mu}{Temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}}},{u_{1} = {\begin{bmatrix} u_{c\; 1} \\ u_{m\; 1} \\ u_{e\; 1} \end{bmatrix}\text{:}\mspace{14mu}{Heater}\mspace{14mu}{electric}\mspace{14mu}{power}}},{{T_{0} = {\begin{bmatrix} t_{0} \\ t_{0} \\ t_{0} \end{bmatrix}\text{:}\mspace{14mu}{Preset}\mspace{14mu}{temperature}\mspace{14mu}{of}\mspace{14mu}{refrigerant}}};}$ and wherein the values of respective elements of constant matrices A₁₁ to A₂₂ and B₁₁ to B₂₂ are determined using the method of least squares.
 14. The control method of claim 13, wherein the Luenberger states observer equations for the first simultaneous linear differential equation include: $\frac{\mathbb{d}z_{2}}{\mathbb{d}t} = {{P\left( {z_{2} - T_{0}} \right)} + {Q\left( {y_{1} - T_{0}} \right)} + {{PL}\left( {y_{1} - T_{0}} \right)} + {{Ru}_{1}\mspace{14mu}{and}}}$ P = (A₂₂ − LA₁₂), Q = (A₂₁ − LA₁₁), R = (B₂₁ − LB₁₁); wherein   $\;{{z_{2} = {\begin{bmatrix} z_{c\; 2} \\ z_{m\; 2} \\ z_{e\; 2} \end{bmatrix}\text{:}\mspace{20mu}{States}\mspace{14mu}{observer}}},{and}}$ L:  any  matrix  of  3  rows  ×  3  columns; wherein from the values of variable z₂ that are calculated in the Luenberger states observer equations, an estimated temperature of the sample ({tilde over (x)}₂) is determined with the following equations: ${{\overset{\sim}{x}}_{2} = {z_{2} + {L\left( {y_{1} - T_{0}} \right)}}},{and}$ ${{\overset{\sim}{x}}_{2} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 2} \\ {\overset{\sim}{x}}_{m\; 2} \\ {\overset{\sim}{x}}_{e\; 2} \end{bmatrix}\text{:}\mspace{14mu}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}}};$ ${{wherein}\mspace{14mu}{the}\mspace{14mu}{following}\mspace{14mu}{unit}\mspace{14mu}{matrix}\mspace{14mu}{is}\mspace{14mu}{used}\mspace{14mu}{as}\mspace{14mu}{the}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{L:L}} = {\begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}.}$
 15. The control method of claim 14, wherein variation of a sample stage temperature is measured; wherein the second simultaneous linear differential equation includes: ${{{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} y_{1} \\ x_{2} \\ x_{3} \end{bmatrix}} = {{\begin{bmatrix} A_{11} & A_{12} & A_{13} \\ A_{21} & A_{22} & A_{23} \\ A_{31} & A_{32} & A_{33} \end{bmatrix}\begin{bmatrix} {y_{1} - T_{0}} \\ {x_{2} - T_{0}} \\ {x_{3} - T_{0}} \end{bmatrix}} + {\begin{bmatrix} B_{11} \\ B_{21} \\ B_{31} \end{bmatrix}u_{1}}}},{wherein}}\mspace{14mu}$ ${y_{1} = {\begin{bmatrix} y_{c\; 1} \\ y_{m\; 1} \\ y_{e\; 1} \end{bmatrix}\text{:}\mspace{14mu}{Sensor}\mspace{14mu}{temperature}}},{x_{2} = {\begin{bmatrix} x_{c\; 2} \\ x_{m\; 2} \\ x_{e\; 2} \end{bmatrix}\text{:}\mspace{14mu}{Temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}}},{x_{3} = {\begin{bmatrix} x_{c\; 3} \\ x_{m\; 3} \\ x_{e\; 3} \end{bmatrix}\text{:}\mspace{14mu}{Sample}\mspace{14mu}{stage}\mspace{14mu}{temperature}}},{u_{1} = {\begin{bmatrix} u_{c\; 1} \\ u_{m\; 1} \\ u_{e\; 1} \end{bmatrix}\text{:}\mspace{14mu}{Heater}\mspace{14mu}{electric}\mspace{14mu}{power}}},{and}$ $\;{{T_{0} = {\begin{bmatrix} t_{0} \\ t_{0} \\ t_{0} \end{bmatrix}\text{:}\mspace{14mu}{Preset}\mspace{14mu}{temperature}\mspace{14mu}{of}\mspace{14mu}{refrigerant}}};}$ wherein suffixes c, m, and e represent the center, middle, and edge regions respectively: and wherein the value of a constant matrix A₁₁ to A₃₃ and B₁₁ to B₃₁ is determined using the method of least squares.
 16. The control method of claim 15, wherein the Luenberger states observer equations for the second simultaneous linear differential equation include: ${{\frac{\mathbb{d}}{\mathbb{d}t}\begin{bmatrix} z_{2} \\ z_{3} \end{bmatrix}} = {{P\begin{bmatrix} {z_{2} - T_{0}} \\ {z_{3} - T_{0}} \end{bmatrix}} + {Q\left( {y_{1} - T_{0}} \right)} + {P\begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right.} \end{bmatrix}} + {Ru}_{1}}},{P = \begin{bmatrix} {A_{22} - {L_{1}A_{12}}} & {A_{23} - {L_{1}A_{13}}} \\ {A_{32} - {L_{2}A_{12}}} & {A_{33} - {L_{2}A_{13}}} \end{bmatrix}},{Q = \begin{bmatrix} {A_{21} - {L_{1}A_{11}}} \\ {A_{31} - {L_{2}A_{11}}} \end{bmatrix}},{R = \begin{bmatrix} {B_{21} - {L_{1}B_{11}}} \\ {B_{31} - {L_{2}B_{11}}} \end{bmatrix}},{z_{2} = \begin{bmatrix} z_{c\; 2} \\ z_{m\; 2} \\ z_{e\; 2} \end{bmatrix}},{and}$ $\mspace{11mu}{{z_{3} = {\begin{bmatrix} z_{c\; 3} \\ z_{m\; 3} \\ z_{e\; 3} \end{bmatrix}\text{:}\mspace{14mu}{States}\mspace{14mu}{observer}}},L_{1},{{L_{2}\text{:}\mspace{14mu}{any}\mspace{14mu} 3\mspace{14mu}{rows}\mspace{14mu} \times \mspace{11mu} 3\mspace{14mu}{columns}\mspace{14mu}{matrix}};}}$ wherein from the values of variables z₂ and z₃ that are calculated in the Luenberger states observer equations, an estimated temperature of the sample ({tilde over (x)}₂) is determined using the following equation: ${\begin{bmatrix} {\overset{\sim}{x}}_{2} \\ {\overset{\sim}{x}}_{3} \end{bmatrix} = {\begin{bmatrix} z_{2} \\ z_{3} \end{bmatrix} + \begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right)} \end{bmatrix}}},{wherein}$ ${{\overset{\sim}{x}}_{2} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 2} \\ {\overset{\sim}{x}}_{m\; 2} \\ {\overset{\sim}{x}}_{e\; 2} \end{bmatrix}\text{:}\mspace{14mu}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}}},{and}$ $\;{{{\overset{\sim}{x}}_{3} = {\begin{bmatrix} {\overset{\sim}{x}}_{c\; 3} \\ {\overset{\sim}{x}}_{m\; 3} \\ {\overset{\sim}{x}}_{e\; 3} \end{bmatrix}\text{:}\mspace{14mu}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{sample}\mspace{14mu}{stage}\mspace{14mu}{surface}\mspace{14mu}{temperature}}};{and}}$ wherein  the  following  unit  matrices are  used  as  the  value  of  L₁  and  L₂: ${L_{1} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}},\;{{{and}\mspace{14mu} L_{2}} = {\begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}.}}$
 17. The control method of claim 16, wherein when temperature of the sample reaches 30° C. for a predetermined amount of time, the sample is processed while controlling the temperature of the sample, under one of three conditions: (1) a condition under which the plasma heat input in the center region is large and a condition under which the preset temperatures of the center/middle/edge regions are controlled to 40/30/30°, respectively, (2) a condition under which the heat input in the middle region is large and a condition under which the preset temperatures of the center/middle/edge regions are controlled to 30/40/30° C., respectively, and (3) a condition under which the heat input in the edge region is large and a condition under which the preset temperatures of the center/middle/edge regions are controlled to 30/30/40° C., the variation of a temperature of a sample (x₂), and sensor temperature (y₁) is measured; wherein the relationship between the temperature of the sample (x₂), sensor temperature (y₁), estimated value ({tilde over (x)}₂) of the temperature of the sample, estimated value of the sample stage surface temperature ({tilde over (x)}₃), and plasma heat input (u₂) are approximated using the following equation: ${{\frac{\mathbb{d}x_{2}}{\mathbb{d}t} = {{A_{41}\left( {y_{1} - T_{0}} \right)} + {A_{42}\left( {{\overset{\sim}{x}}_{2} - T_{0}} \right)} + {A_{43}\left( {{\overset{\sim}{x}}_{3} - T_{0}} \right)} + {A_{44}\left( {x_{2} - T_{0}} \right)} + {B_{42}u_{2}}}},{wherein}}\mspace{11mu}$ $\;{{u_{2} = {\begin{bmatrix} u_{c\; 2} \\ u_{m\; 2} \\ u_{e\; 2} \end{bmatrix}\text{:}\mspace{20mu}{Amount}\mspace{14mu}{of}\mspace{14mu}{plasma}\mspace{14mu}{heat}\mspace{14mu}{input}}};}$ wherein a constant matrix A₄₁ to A₄₄ and B₄₁ is determined using the method of least squares; and wherein it is not necessary to accurately measure the amount of plasma heat input u₂, because under each condition, the respective one of the following three types of independent matrices is used: ${{u_{2}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu} 1} = \begin{bmatrix} 1 \\ 0 \\ 0 \end{bmatrix}},{{u_{2}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu} 2} = \begin{bmatrix} 0 \\ 1 \\ 0 \end{bmatrix}},{{u_{2}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu} 3} = {\begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix}.}}$
 18. The control method of claim 17, wherein after the amount of plasma heat input u₂ under an actual process condition is calculated, the estimated value ({tilde over (x)}₂) of the temperature of the sample in plasma heat input is measured using the following Luenberger states observer equations: ${\frac{\mathbb{d}z_{4}}{\mathbb{d}t} = {{P^{\prime}\begin{bmatrix} {z_{2} - T_{0}} \\ {z_{3} - T_{0}} \\ {z_{4} - T_{0}} \end{bmatrix}} + {Q^{\prime}\left( {y_{1} - T_{0}} \right)} + {P^{\prime}\begin{bmatrix} {L_{1}\left( {y_{1} - T_{0}} \right)} \\ {L_{2}\left( {y_{1} - T_{0}} \right)} \\ {L_{3}\left( {y_{1} - T_{0}} \right)} \end{bmatrix}} + {R^{\prime}\begin{bmatrix} u_{1} \\ u_{2} \end{bmatrix}}}},{and}$ ${P^{\prime} = \left\lbrack {A_{42} - {L_{3}A_{12}A_{43}} - {L_{3}A_{13}A_{44}}} \right\rbrack},{Q^{\prime} = {A_{41} - {L_{3}A_{11}}}},{R^{\prime} = \left\lbrack {{- L_{3}}B_{11}B_{42}} \right\rbrack},{z_{4} = {\begin{bmatrix} z_{c\; 4} \\ z_{m\; 4} \\ z_{e\; 4} \end{bmatrix}\text{:}\mspace{14mu}{States}\mspace{14mu}{observer}}},{{L_{3}\text{:}\mspace{14mu}{any}\mspace{14mu}{matrix}\mspace{14mu}{of}\mspace{14mu} 3\mspace{14mu}{rows}\mspace{14mu} \times 3\mspace{14mu}{columns}};{and}}$ x̂₂ = z₄ + L₃y₁, wherein ${\hat{x}}_{2} = {\begin{bmatrix} {\hat{x}}_{c\; 2} \\ {\hat{x}}_{m\; 2} \\ {\hat{x}}_{e\; 2} \end{bmatrix}\mspace{14mu}{Estimated}\mspace{14mu}{value}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}{in}\mspace{14mu}{the}}$ presence  of  the  plasma  heat  input; wherein the value of the amount of plasma heat input u₂ is changed from zero to a value obtained simultaneously at the start of plasma heat input; wherein the following unit matrix is used as a value of L₃: $L_{3} = {\begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}.}$
 19. The control method of claim 18, wherein a transformation matrix from sensor temperature to temperature of the sample is determined by using the following equation, using the amount of change in temperature of the sample and sensor temperature: ${{{{{C = {{\begin{bmatrix} {\Delta\; x_{c\; 1}} & {\Delta\; x_{c\; 2}} & {\Delta\; x_{c\; 3}} \\ {\Delta\; x_{m\; 1}} & {\Delta\; x_{m\; 2}} & {\Delta\; x_{m\; 3}} \\ {\Delta\; x_{e\; 1}} & {\Delta\; x_{e\; 2}} & {\Delta\; x_{e\; 3}} \end{bmatrix}\begin{bmatrix} {\Delta\; y_{c\; 1}} & {\Delta\; y_{c\; 2}} & {\Delta\; y_{c\; 3}} \\ {\Delta\; y_{m\; 1}} & {\Delta\; y_{m\; 2}} & {\Delta\; y_{m\; 3}} \\ {\Delta\; y_{e\; 1}} & {\Delta\; y_{e\; 2}} & {\Delta\; y_{e\; 3}} \end{bmatrix}}^{- 1}{{wherein}\begin{bmatrix} {\Delta\; x_{c\; 1}} \\ {\Delta\; x_{m\; 1}} \\ {\Delta\; x_{e\; 1}} \end{bmatrix}}{{and}\begin{bmatrix} {\Delta\; y_{c\; 1}} \\ {\Delta\; y_{m\; 1}} \\ {\Delta\; y_{e\; 1}} \end{bmatrix}}{are}\mspace{14mu}{amounts}\mspace{14mu}{of}\mspace{14mu}{change}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu}(1)\mspace{14mu}{of}}}{temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}\mspace{14mu}{and}\mspace{14mu}{sensor}\mspace{14mu}{temperature}\mspace{14mu}{respectively}},{\begin{bmatrix} {\Delta\; x_{c\; 2}} \\ {\Delta\; x_{m\; 2}} \\ {\Delta\; x_{e\; 2}} \end{bmatrix}{{and}\begin{bmatrix} {\Delta\; y_{c\; 2}} \\ {\Delta\; y_{m\; 2}} \\ {\Delta\; y_{e\; 2}} \end{bmatrix}}{are}\mspace{14mu}{amounts}\mspace{14mu}{of}\mspace{14mu}{change}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu}(2)\mspace{14mu}{of}}}{temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}\mspace{14mu}{and}\mspace{14mu}{sensor}\mspace{14mu}{temperature}\mspace{14mu}{respectively}},{\begin{bmatrix} {\Delta\; x_{c\; 3}} \\ {\Delta\; x_{m\; 3}} \\ {\Delta\; x_{e\; 3}} \end{bmatrix}{{and}\begin{bmatrix} {\Delta\; y_{c\; 3}} \\ {\Delta\; y_{m\; 3}} \\ {\Delta\; y_{e\; 3}} \end{bmatrix}}{are}\mspace{14mu}{amounts}\mspace{14mu}{of}\mspace{14mu}{change}\mspace{14mu}{under}\mspace{14mu}{condition}\mspace{14mu}(3)\mspace{14mu}{of}}}{temperature}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu}{sample}\mspace{14mu}{and}\mspace{14mu}{sensor}\mspace{14mu}{temperature}\mspace{14mu}{{respectively}.}$ temperature of the sample and sensor temperature respectively.
 20. The control method of claim 18, wherein the increase in temperature of the sample (Δx) caused by plasma heat input is calculated from the increase in sensor temperature (Δy), by using the following equation: Δ x = C ⋅ Δ y  wherein ${{{Increase}\mspace{14mu}{of}\mspace{14mu}{wafer}\mspace{14mu}{temperature}\mspace{14mu}\Delta\; x} = \begin{bmatrix} {\Delta\; x_{c\;}} \\ {\Delta\; x_{m\;}} \\ {\Delta\; x_{e}} \end{bmatrix}},{{{Increase}\mspace{14mu}{of}\mspace{14mu}{sensor}\mspace{14mu}{temperature}\mspace{14mu}\Delta\; y} = {\begin{bmatrix} {\Delta\; y_{c\;}} \\ {\Delta\; y_{m}} \\ {\Delta\; y_{e}} \end{bmatrix}.}}$ 